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Abstract. In a recent article the first three authors proved that in dimension 4m + 1 all homo- 
topy spheres that bound parallelizable manifolds admit Einstein metrics of positive scalar curvature 
which, in fact, are Sasakian-Einstein. They also conjectured that all such homotopy spheres in di- 
mension Am — 1, m > 2 admit Sasakian-Einstein metrics [BGK03], and proved this for the simplest 
case, namely dimension 7. In this paper we describe computer programs that show that this conjec- 
ture is also true for 11-spheres and 15-spheres. Moreover, a program is given that determines the 
partition of the 8610 deformation classes of Sasakian-Einstein metrics into the 28 distinct oriented 
diffomorphism types in dimension 7. 



1. Introduction 

In a recent article the first three authors gave a method for constructing Einstein metrics of positive 
scalar curvature on odd-dimensional homotopy spheres [BGK03]. By Kervaire and Milnor [KM63] 
and Smale [Sma62], for each n > 5, differentiable homotopy spheres of dimension n form an Abelian 
group 0„, where the group operation is connected sum. 0„ has a subgroup bP n+ i consisting of those 
homotopy n-spheres which bound parallelizable manifolds V n +i- Kervaire and Milnor [KM63] proved 
that W-2m+i = for m > 1, bP^ m+ 2 = 0, or Z2 and is Z2 if 4m + 2 ^ 2* — 2 for any i > 3. The most 
interesting groups are bP^ m for m > 2. These are cyclic of order 



bP im \ = 2 Im - I (2' lm - 1 - 1) numerator 



4P„ 



where B m is the m-th Bernoulli number. Thus, for example |6P 8 | = 28, |&P 12 | = 992, |6P 16 | = 8128 
and |6P 20 | = 130,816. In the first two cases these include all exotic spheres. The correspondence is 
given by 

KM : S » |r(l/ 4m (S)) mod \bP 4m \, 

where V4 m (S) is any parallelizable manifold bounding E and t is its signature. Let Ej denote the 
exotic sphere with 1£"M (Ej) = i. 

In [BGK03] the authors proposed the following: 

Conjecture 1. The construction of [BGK03] yields Einstein metrics on every exotic sphere that 
bounds a parallelizable manifold. 

The construction is described in sections 2-3. The method gives Einstein metrics whose isometry 
group is one dimensional and they are even Sasakian-Einstein. 

In [BGK03] the conjecture was shown to be true in dimensions Am + 1. In dimension 7 we were 
also able to verify it; the relevant signature calculations were carried out by a computer. 

The main aim of this paper is to provide more evidence for our conjecture by demonstrating that 
it is true in dimensions 11 and 15 as well. More precisely we show 

Theorem 2. Every homotopy sphere Ej £ bPn and E, G 6P16 admits at least one Einstein metric. 



We also give a complete enumeration of all oriented diffeomorphism types in dimension 7, namely, 
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Theorem 3. In dimension 7, £j admits at least rn inequivalent deformation classes of Einstein met- 
rics, where (ni,...,n 28 ) = (376,336,260,294,231,284,322,402,317,309,252,304,258,390,409,352, 
226, 260, 243, 309, 292, 452, 307, 298, 230, 307, 264, 353), giving a total of 8610 cases. 

Actually for dimensions 11 and 15, just as in dimension 7, we do get several deformation types, 
but the signature was computed only for a sample of all cases. For instance, in dimension 15 our 
method gives at least 10 50 deformation classes of Einstein metrics on all homotopy 15-spheres, and 
even their complete enumeration is impossible with the current programs and facilities. 

2. Brieskorn-Pham Singularities and Their Links 

For a = (ai, . . . , a m ) e set F a (z) := E"=i z i i ■ Consider a Brieskorn-Pham singularity 

Y(a) := (F a (z) = 0) C C m , and its link L(a) := Y(a) n S 2m_1 (l). 

Set C — lcm(a i : i = 1, . . . , m). Both Y(a) and L(a) are invariant under the C*-action 

{z u ...,z m ) i ► (A c/ai z 1 ,...,A c / a '"z m ). 

If we denote w = (wi, . . . ,w m ) = (C/ai, . . . ,C/a m ), then F a is a weighted homogeneous polynomial 
on C m with weight w and degree C, i.e., 

F gL (X Wl zi, X Wrn z m ) = X c F 8L (zi, ...,z m ). 

Consider the orbit spaces: X orb (a) := Y(a) \ {0}/C* and the weighted projective space P(w) := 
(C m \ {0}/C*). We get a commutative diagram 

L(a) > S 2 ™- 1 

X orb (a) > P(w). 

It is known that the sphere S* 2 ™ -1 can be given a Sasakian structure with respect to the projection 

g2m-i >p( w ) associated to the characteristic foliation [YK84]. In such a case the embedding 

L(a) ^S* 2 " 1 " 1 is Sasakian and X orb (a) is the horizontal space of the characteristic foliation of the 

link L(a) [BG01]. 

3. Orbifolds and Einstein Metrics 

Let C- 7 = lcm(ai : i ^ j), bj = gcd(cij , C 3 ) and dj = aj/bj. The following result was estab- 
lished in [BGK03]. 

Theorem 4. The orbifold X orb (a) — Y(a)/C* is Fano and has a Kdhler-Einstein metric if 

(1) i < rr=i h 

(2) E™ii<l + ^min,{^-} J and 

(3) Er=i^<l + ^niin i , j { 5 l-}. 

In this case the link L(a) admits a Sasakian- Einstein metric with one- dimensional isometry group. 

The first inequality is necessary for X orb (a) to be Fano. Hence, it is also necessary for the link 
L(a) to admit any Sasakian-Einstein structure. The second inequality is necessary for our algebraic 
approach to Kahlcr-Einstein metrics to work, while the third inequality is most likely an artifice of 
our estimates. Hopefully, it is not needed at all. We should reiterate that the failure of our method 
does not imply that X orb (a) cannot admit a positive Kahler-Einstein metric as long as X orb (a) is 
Fano. 
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For any m > 3 there are infinitely many m-tuples satisfying the conditions of Theorem 4. For 
example, we can take a = (m — 1, . . . , m — 1, k), where gcd(m — 1, k) = 1 and k > (m — l)(m — 2). 
However, in this paper we are only interested in the case when the link L(a) is a homotopy sphere 
and, as we shall see, L(m — 1, . . . , m — 1, k) is not. 

4. Homotopy Spheres as Brieskorn-Pham Links 

To every m-tuple a one can associate a graph G(a) whose m vertices are labeled by d\, ■ ■ ■ , a m . 
Two vertices dj and a, are connected if and only if gcd(aj,Oj) > 1. Let C ev denote the connected 
component of G(a) determined by the even integers. Note that all even vertices belong to C ev , but 
C ev may contain odd vertices as well. Brieskorn shows that: 

Theorem 5. [Bri66] The link L(a) (with m > 3) is a homotopy sphere if and only if either of the 
following hold: 

(1) G(a) contains at least two isolated points, or 

(2) G(a) contains one odd isolated point and C ev has an odd number of vertices and for any 
distinct di , CLj G Cev, gcd(ai,a,j) = 2. 

We observe that, in each dimension, there are only finitely many m-tuples that yield homotopy 
spheres and satisfy the conditions (4.1-3). For that we introduce 

Example 6 (Euclid's or Sylvester's sequence). (See [GKP89, Sec.4.3] or [Slo03, A000058].) 
Consider the sequence defined by the recursion relation 

Cfe+i = ci • • • c fc + 1 = c\ - c k + 1 

beginning with c\ = 2. We call this sequence the extremal sequence. It starts as 

2, 3, 7, 43, 1807, 3263443, 10650056950807, 

and it is easy to see (cf. [GKP89, 4.17]) that 

"1 1 1 
2J- = 1 — = 1 . 

c i c m+l J- Ci • • • C m 

In [Sou03] it was proved that if the sum of reciprocals of m natural numbers is less than 1, then it is 
at most 1 — l/(c m+ i — 1). Thus, in this sense the sequence {ci} is extremal. 

We use the sequence Cj to show that the number of m-tuples that yield homotopy spheres and 
satisfy the conditions (4.1-3) is finite. Without loss of generality we shall assume that the exponents 
are arranged in the non-decreasing order. 

Proposition 7. Assume that a € Z™ satisfies the conditions (4-1-2) and (5). Then au < (m — k + 
l)(cfe — 1), for k = 1, . . . ,ra — 1 and a m < -^z^icm — !)• In particular, the number of such m-tuples 
is finite for each m > 3. 

Proof. Step 1. We first observe that Y^i=\ — < 1- For otherwise we would have 

1 1 m- 1 1 2 
+ — < • — < — , 

which is impossible. 

Step 2. Now, assume that ^- < 1. Then it is also < 1 — l/(ck+i — 1). The remaining m — k 

reciprocals must sum to more than l/(cfe + i — 1), hence we obtain that ak+i < (m — k)(ck+i — 1). By 
Step 1 this takes care of all a% for i < m — 1 and also of a m if X^i* ~ < !■ 

Thus we are left with Step 3: ^"ii 1 > 1 . If equality holds there is no bound for a m ; however, 
in this case £(a) is not a homotopy sphere, since Theorem 5 says that at least one of the a\, . . . , a m _\ 
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(or half of it) is relatively prime to the others, and this implies that we cannot get an integer as a 
sum of reciprocals. Otherwise we have 

-l 



E 1 



> 1 + > 1 + — ; = 1 



r[ a% a\---a m -\ ml ■ ci ■ ■ ■ c m -! ml ■ (c m - 1) ' 

Thus we obtain that 



1 m 1 

— <E 1 

n t—i n . 



1 1^1 to- 1 1 

< 1 



ml ■ (c m — 1) a m ■f—' a; m — 2 a 



Comparing the two ends gives that 



a m < ^(cm - 1). 

m — 2 



□ 



We wrote a simple program which we call candidates . c . This is a C code which enumerates all 
ordered m-tuples satisfying the conditions (4.1-3) and (5.1) or (5.2) in any given range af ln < m < 
a™, i — 1, . . . , m, with the condition that a™ m < • • < a™. In principle, for any m > 4, the 
program can be used to enumerate all m-tuples of this type. However, this is not feasible already 
for m = 7. On the other hand, the program has the flexibility to "hunt" for such m-tuples in any 
specified region of the integral lattice defined by Proposition 7. 

5. DlFFEOMORPHISM TYPES - BRIESKORN, ZAGIER, AND HlRZEBRUCH 

By Theorem 5, we know when L(a) is a homotopy sphere. We now would like to be able to 
determine the diffeomorphism types of various links. In this article, we are only interested in the 
case when m = 2k + 1. 

In this case, the diffeomorphism type of a homotopy sphere £(a) £ W"2m-2 is determined [KM63] 
by the signature r(M) of a parallelizable manifold M whose boundary is £^ m ~ 3 . By the Milnor 
Fibration Theorem [Mil68] we can take M to be the Milnor fiber M^" l ~ 2 which, for links of isolated 
singularities coming from weighted homogeneous polynomials is diffeomorphic to the hypersurface 
{z e C m | F(zi, ■■■ ,z m ) = 1}. 

Brieskorn shows that the signature of M 2m_2 can be written combinatorially as 

2k 

T{Mt k ) = #{x G Z 2fe+1 I < x, < a t and < ^ — < 1 mod 2} 

- n a i 

5.1 

Ik 

- #{x G Z 2fc+1 | < Xi < a t and 1< ^ — < 2 mod 2}, 
where m = 2k + 1. 

Using a formula of Eisenstein, Zagier (cf. [Hir71]) has rewritten this formula as: 

5.2 ,(Mf ) = til V cot cot ... cot + 

N 2N 2 ao 2a 2k 

where N is any common multiple of the a^'s. Both formulas are quite well suited to computer use. 
We wrote a second C code which we call sig.c. For any m-tuple with m = 2k + 1 = 5,7, 9 sig.c 
computes the signature r(a) := r(M^ fc ) and the diffeomorphism type of the link using either of the 
above formulas. Furthermore, one can use sig.c to compute signature and diffeomorphism type of 

lr The codes candidates. c and sig.c as well as all the relevant data files mentioned later can be downloaded at 
http://www.math.unm. edu/"galicki/papcrs/codcs. html. 
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a single ra-tuple, or one can select an arbitrary set of m-tuples X and compute the signature and 
diffcomorphism type associated to every m-tuple a El. One last feature of sig. c is that, provided 
an appropriate option is chosen, the program will start computing diffcomorphism type g(a) of each 
m-tuple a G X until it finds all possible oriented diffcomorphism types in bP2 m -2 after which it stops. 



6. The proofs 

Proof. [Theorem 3] In dimension 7 the candidates . c can be run in the maximal range specified 
by Proposition 7. The result is exactly 8610 solutions. These solutions become an input data file X 
for the signature computation using sig. c with either Bricskorn or Zagicr formula. In the case of 5- 
tuplcs the choice is not important. The signature computation takes a couple of hours on a Pentium 4 
processor and the result is a list of 8610 5-tuplcs a = (a 1; a 2 , a 3 , a 4 , a 5 ) each with a number g(a) E Z 2 s 
which determines the oriented diffcomorphism type of L(a). The results are contained in the output 
file 7spheres.txt. This file can be easily sorted grouping 5-tuples with the same g(a) and we get 
the result described in Theorem 3. □ 

Proof. [Theorem 2] In dimension 11 the candidates. c cannot be run in the maximal range of 
Proposition 7. The complete enumeration would take too long a time. Instead, the code candidates . c 
is used to select 7-tuples in a specified range. This will become an input file X for the subsequent 
signature computation. One important point in selecting X is that C = lcm(ai, . . . , a-j) should not be 
too large. The time of every individual signature computation with sig. c is approximately linear in 
C. Another relevant point is that bP\i = Z 99 2 so that \X\ should be sufficiently large. For example, 
we can ask candidates . c to search for 7-tuples in the following range: 2 < a\ < 6, 3 < a 2 < 11 and 
i + 1 < ai < 30 for i = 3, 4, 5, 6, 7. This guarantees a relatively small C < 66 • 30 5 for all solutions 
and \X\ = 21, 535. One should point out that there is nothing special about the choice of X - other 
choices can be equally successful in yielding the desired result. We now want to determine if we 
find all g(a) E Z992 among a G X. This is done by feeding each 7-tuple a G X into sig.c with the 
following option: the program will calculate the signature r(a) and the diffeomorphism type g(a) of 
each 7-tuple a £ I in the order specified by X. Any 7-tuple ael with a diffcomorphism type g(a) 
not previously found gets automatically recorded into the output file A. Once the program finds all 
992 oriented diffcomorphism types it stops. The output file contains a subset of the original input file 
(hopefully) containing exactly 992 7-tuples. All this work can done on a single PC with a Pentium 4 
processor. An example of an output file A called llspheres . txt can be found at the URL mentioned 
in the earlier footnote. We needed approximately 9000 7-tuples to find the 992 necessary to prove 
Theorem 2. 

In dimension 15 we repeat the steps outlined in the 1 1-dimcnsional case. Selecting appropriately 
large data file with candidates, c is not a problem. This can be done on a single PC. Given that 
bPi & ~ Z 8 i28 one needs an input file X with about 80K 9-tuplcs for the signature computation with 
sig. c. More challenging problem has to do with computing signature for these 9-tuplcs. To minimize 
computing time some care should be given to how X is selected. The length of a single computation 
varies depending on (1) a = (ai, . . . , ag) itself; (2) the formula used for the signature computation; 
(3) the processor's speed. Also, this is an easy parallelization task because it consists of tens of 
thousand runs which are almost completely independent of each other. The only co-ordination that 
is required is to stop the process when all g(a) G Z 8 i 2 8 are found. 

We actually generated two sets X\ , X2 for the signature calculation. The first set X\ was was created 
by appropriately restricting the size of all exponents. After the calculations for X\ were completed 
a second set X 2 was chosen to select 9-tuples with a restricted upper bound on C — lcm(ai, . . . , ag). 
We first used the Zagier formula 5.2 to calculate the signature r(a) of the 9-tuples in the selected 
input files Xi 7 X 2 . Zagier's formula was chosen as this calculation is much faster for most individual 
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a's. Exactly how much faster depends on the ratio a\ ■■ -ag/C. If C = a\ ■ --ag then the Bricskorn 
formual 5.1 is slightly faster. On the other hand, the problem with using the Zagier formula for very 
large C is that there is a large round-off error on Intel x86 processors even at maximum precision. 
When C is of the order of 10 9 this error becomes large enough that g(a) is sometimes calculated 
incorrectly. While this was not an issue for all 5-tuples and also for carefully selected 7-tuples the 
case of 15-spheres was more of a problem. Instead of forcing the program to do a better round-off 
error control with the Zagier option, we decided to do the first calculation with the Zagier formula 
and then verify all signature calculations for the candidate solution with the Bricskorn formula 5.1. 
By its nature, this formula does not have any round-off error. At the end we actually generated 
two disjoint sets of 9-tuples. One is contained in the file 15spheresA.txt. The other one is in 
15spheresB.txt. □ 

The Zagier calculation on the first set Ti was done at the University of Melbourne on an IBM eS- 
erver 1350 which is a cluster of 48 2.4GHz Intel Xeon processors. The calculation leading to the data 
set 15spheresA . txt took approximately 9500 hours of processor time and tested nearly 70,000 candi- 
dates. The Bricskorn verification was performed on 15spheresA . txt at the University of New Mexico 
High Performance Computer Center on a 256 node cluster of 733Mhz processors. This required 80,000 
hours of processor time. In the case of one 9-tuple the code calculating with the Zagier formula yielded 
the wrong answer: g z (3, 4, 8, 8, 9, 43, 83, 85, 97) = 3323 while g B (Z, 4, 8, 8, 9, 43, 83, 85, 97) = 3322 is 
correct. Note that for this particular example C a = 2,118,701,160. It is in the 10 9 range where 
sig.c becomes unreliable with the Zagier option. An additional search for a 9-tuplc with that par- 
ticular oriented diffeomorphism type was performed so that 15spheresA.txt actually contains the 
full set of 8128 examples. We replaced it with a = (6, 6, 6, 6, 6, 10, 25, 59, 73) which came out of T^. 
Note that here C a = 646, 050 which is smaller by 3 orders of magnitude. 

Realizing that one can do much better by a careful selection of candidates with low C a = 
lcm(ai, . . . , ag) we used candidates . c to select more "efficient" input data set Ti- As a result it was 
possible to obtain all 8128 distinct 17(a) 's calculating with the Zagier option in only about 160 hours 
on the University of Melbourne facility. A very significant improvement indeed. That second calcu- 
lation generated 15spheresB.txt. The Brieskorn verification was performed on 15spheresB.txt at 
the University of Melbourne facility and it took only 1700 hours. No errors were found in the Zagier 
calculation which is no surprise: a typical C a for 9-tuples of I2 was about 3 orders of magnitude 
lower. 

Note that one can easily improve the "least one" statement of Theorem 2 by repeating the same 
calculation with several disjoint input files T. It is a simple exercise to do it for 7-tuples and much 
more time consuming in the case of 9-tuples. For 9-tuples we actually showed that there are at least 
two Sasakian-Einstein metrics on each homotopy sphere a} 5 G bPie as the lists 15spheresA.txt and 
15spheresB.txt arc disjoint. 

On the other hand, to calculate signatures of all candidate 7-tuples and 9-tuples to get the state- 
ment similar to the one expressed in Theorem 3 would take thousands of years with the present 
technology. 

Remark 8. It is clear that our approach breaks down for (2n+ l)-tuples, where n is "large enough" . 
What is exactly "large enough" depends on several factors. Our rough estimate indicates that 
assuming the same facilities and the same codes are used it would take about 100 years to do the 
same calculation for 19-spheres. No doubt the sig.c code can be improved to calculate faster. On 
the other hand an average C for 11-tuples will be at least 10 3 larger that in the 9-tuple case. In 
addition, bP 2 o = 130, 816 is much bigger. Taking these two factors into account, a calculation for 
19-spheres would take about 10 4 times longer that a similar calculation for 15-spheres. 
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